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ABSTRACT 



An error probability analysis is performed for a conventional noncoherent M-ary 
orthogonal frequency-shift keying (MFSK) receiver employing fast frequency-hopped 
(FFH) spread spectrum waveforms transmitted over a frequency-nonselective, slowly 
fading Ricean channel with partial-band noise interference . Each diversity reception is 
assumed to fade independently. The partial-band interference is modeled as a Gaussian 
process. The effects of wideband thermal noise are also included. The energy per hop is 
held constant; thus, as diversity increases, energy per symbol increases. Previous 
analyses considered only constant energy per symbol systems, however, practical military 
systems are likely to employ fixed hop rates. There is some performance enhancement to 
be obtained from implementing diversity in a conventional FFH/MFSK system with fixed 
hop rates, but partial band interference still results in significant degradation. 

Additionally, the performance of this FFH receiver is investigated over the seime 
channel in the presence of partial-band tone jamming without diversity for the case of 
binary frequency-shift keying (BFSK) when both the signal and the jammer can fade 
independently. Performance when only a single jamming tone per hop slot is allowed is 
compared to that obtained when two jamming tones per hop slot are possible. When the 
jamming signal experiences Rayleigh fading there is very little degradation of the 
jammer's effectiveness as compared to when the jamming signal is not affected by fading. 
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I. INTRODUCTION 



A fast frequency-hopping (FFH) communication system is a subset of spread 
spectrum communications that utilizes a bandwidth greatly exceeding that required for 
the information signal alone. Frequency-hopping spread spectrum is fundamentally 
different from direct sequence (DS) spread spectrum in the technique of signal generation 
and recovery. This thesis focuses on frequency-hopping spread spectrum systems 
because of their practical military importance. As spread spectrum systems grow more 
popular and occupy wider communications bandwidth, the likelihood of both hostile and 
non-hostile sources of narrowband interference also grows. It is important for mission 
planners to have a reliable estimate of the degradation that their commimication systems 
will suffer under partial-band noise interference as well as tonal interference. Fast 
frequency-hopping spread spectrum techniques have evolved to counter the threat of 
intentional jamming [Ref 1]. This thesis presents an error probability analysis of the 
conventional fast frequency-hopped orthogonal M-ary frequency-shift keying 
(FFH/MFSK) receiver with noncoherent detection for communications over channels 
with Ricean fading of the signal and partial-band noise interference. The influence of 
partial-band tonal interference on the binary orthogonal frequency-shift keying 
(FFH/BFSK) receiver without diversity where both the signal and jammer may fade 
independently is also considered. 
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A. FFH/MFSK 



Frequency-shift keying is popular as a signaling scheme because it allows for 
noncoherent reception of the signal. A typical MFSK signal set can be expressed as 



5 1 (t) = Jl A ccos(2nfot -I- 0 1 ) 

52(0= j2AcCOs[2n(fo + Af)t + B2] 

and so on to 

sm(0= J 2 AcCos[2n(fo + (M- l)A/)r-i-0M] 



( 1 ) 

( 2 ) 

(3) 



whereto is the lowest signal tone, A/is the tone spacing, and Qm is the phase associated 
with each tone. This set is then modulated by a carrier that varies pseudorandomly in 

frequency for transmission as a frequency-hopped signal. 

All communication signals suffer from interference and noise. One method to 
overcome this degradation is to transmit the signal more than once thus providing a form 
of diversity. Fast frequency-hopping employs this redundancy as well as deliberately 
spreading the bandwidth. Frequency-hopping systems in which several symbols are 
transmitted per hop are considered slow frequency-hopping, while those transmitting 
several hops per symbol are considered fast frequency-hopping. Our FFH/MFSK 
transmitter performs L hops per data symbol, which results in a diversity of L. At the 
receiver the dehopped signals are recovered noncoherently by two correlators in phase 
quadrature with the dehopped signal waveform. The correlator outputs are sampled every 
Tf, seconds, where is the hop period. The sampled output of each correlator pair is 
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squared and then these outputs are summed L times to obtain the decision statistic for 
each branch of the M-ary detector. The largest signal detected is selected as the 
transmitted symbol. A typical receiver structure is shown in Figure 1. 

Proper reception and demodulation of the spread signal depend on the recovery of 
several pieces of timing information. Both the sender and receiver need the same 
pseudorandom sequence operating synchronously. Also required are the symbol period 
and the hop period. In practice these are estimated from the received waveforms. This 
thesis assumes that this information is recovered without error. 

B. MULTI-PATH EFFECTS 

The losses experienced by the signal during propagation is worthy of an entire study 
in itself However, it is useful to make some general observations on the composition of 
the received signal. It is possible and even likely that the received signal arrives at the 
receiver after transiting a variety of different paths. Signals traveling a longer distance 
arrive delayed relative to the direct path signal. This leads to multi-path effects. The 
magnitude of the multi-path effects depends on the magnitude of the delayed signal 
strength versus the direct path signal strength. It is common to consider the sum of all 
delayed signals as the diffuse component of the received signal. [Ref 2] 
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Figure 1. Typical MFSK Receiver Structure 
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Choose Largest 




The consequence of multi-path reception is to cause the signal to fade in a time 
varying fashion. We can broadly characterize the channel conditions by examining the 
magnitude of the direct signal power to the diffuse signal power. The Ricean channel is 
the general case. Channels that have nearly all the received signal energy in the direct 
component, i.e., direct-to-diffuse ratio greater than 100, have essentially no fading. In the 
limit, an infinite direct-to-diffuse ratio implies a Gaussian channel. Channels that have 
nearly all the received signal energy in the diffuse component, i.e., direct-to-diffuse ratio 
less than one, have strong fading. A direct-to-diffuse ratio of zero implies a Rayleigh 
channel. For direct-to-diffuse ratios between these extremes the channel experiences 
Ricean fading. This thesis examines the performance in each of these broad categories. 

In addition to these broad categories, the time varying nature of the channel can be 
described as slow or fast. In this thesis the channel properties are assumed to be constant 
over the duration of a hop and, therefore, slowly varying. Further, the channel may 
introduce some signal distortion arising from the treatment of sinusoids comprising the 
signal set within a hop differently. This distortion is characteristic of frequency selective 
channels. However, it is reasonable to assume that the signal sinusoids experience the 
stune multi-path effects. This is the case in frequency nonselective channels. One 
measure of this phenomenon is the coherence bandwidth of the signal. The coherence 
bandwidth is the frequency range over which the signal frequencies pass through without 
distortion. This can be summarized mathematically as 

(A/), = (4) 

^ m 
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where (A/)^ is the coherence bandwidth and T„, is the multi-path spread of the channel. 

A frequency nonselective cannel displays a coherence bandwidth that is larger than the 
signal bandwidth. The rate of fading is related to 

(A/). = ^ (5) 

The coherence time is (At);, and is the Doppler spread of the channel. Slowly 
fading chaimels display a large coherence time or, conversely, a small Doppler spread. 
These descriptors are discussed further in the system description. [Ref. 2] The use of 
diversity to mitigate the multi-path effects for conventional MFSK has been widely 
investigated [Refs. 2, 3, 4]. 

When analyzing performance it is important to distinguish between a constant energy 
per hop system and a constant energy per symbol system. As diversity increases, the total 
symbol energy in a constant energy per hop system increases, while the symbol rate 
decreases. However, in the constant energy per symbol system, increasing diversity, L, 
implies decreasing hop duration, T,,. Hence, a constant data rate is maintained, but the 
energy per hop is reduced. Since many practical military communication systems employ 
a fixed hop rate and a variable data rate, the constant energy per hop assumption is more 
logical [Ref. 5]. 

C. CONSTANT ENERGY PER HOP SYSTEMS 

Previous work has examined the performance of the noncoherent MFSK receiver in a 
Ricean fading channel with partial-band noise interference and constant energy per 
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symbol [Ref. 3]. Additionally, the performance of the conventional BFSK receiver under 
partial-band noise jammin without the multi-path effect of fading is analyzed in [Ref. 6]. 
However, these investigations do not consider constant energy per hop signaling. The 
performance of several different types of diversity combining receivers, including the 
conventional receiver, utilizing constant energy per hop is simulated for Rayleigh fading 
channels in [Ref. 7]. 

Motivating this thesis is the uncertain degree of improved performance offered by 
more elaborate receiver designs, such as the noise normalized receiver, over the 
conventional receiver for variable data rate systems. [Refs. 3, 4] The expense and 
complexity of a more elaborate receiver may not be justified in some circumstances when 
utilizing a constant energy per hop system. 

The constant energy per hop scheme is of practical value because it allows the 
potential for an adaptive signaling scheme in which the sender and receiver can optimize 
the data transmission rate. In poor environments, diversity can be increased at the 
expense of lowered data rate. In favorable environments, the level of diversity, L , can be 
lowered to accommodate a higher data rate. These adaptations will ideally not require 
any hardware modifications and will be transparent to a channel observer. With this 
adaptive scheme in mind , this thesis examines the improvement offered by varying the 
level of diversity possible when the jammer is sub-optimal. 
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D. TONE INTERFERENCE 



Another type of narrowband interference is tone interference. This can consist of a 
single interfering or multiple interfering tones. In this thesis, a performance analysis for a 
FFH/BFSK receiver without diversity over fading channel conditions similar to those 
assumed for noise interference is considered. Since the interfering tone or tones are 
signal-like in nature, they too can suffer multi-path effects. This analysis considers the 
effects of fading on both the signal and jammer. In a FFH/BFSK system an intelligent 
jammer can potentially cause more degradation by splitting his available power over both 
the signal tones. The degree of communication impairment of the single tone interference 
per hop versus two interference tones per hop strategy is also considered. Clearly, the 
greatest performance degradation occurs when the interfering tones correspond exactly to 
the various frequency-hopped symbol tones. Tone jamming where the tones do not 
correspond exactly to the various frequency-hopped symbol tones are not considered in 
this thesis. 
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II. PARTIAL-BAND NOISE INTERFERENCE 



A. SYSTEM DESCRIPTION 

The partial-band noise interference, either intentional or unintentional, considered in 
this thesis is modeled as additive Gaussian noise and, when present, is assumed to be in 
all branches of the MFSK demodulator for any reception of the dehopped signal. 

Thermal noise and other wideband interferences which are also assumed to corrupt the 
signal are modeled as additive white Gaussian noise. Only the signal is assumed to be 
affected by channel fading. The smallest spacing between frequency hop slots is assumed 
larger than the coherence bandwidth of the channel, hence, each dehopped signal fades 
independently [Refs. 2,8]. As discussed in the Introduction, the signal bandwidth is 
assumed to be much smaller than the channel coherence bandwidth, and the channel 
coherence time is assumed to be much larger than the hop duration or, equivalently, the 
hop rate is assumed to be large compared to the Doppler spread of the channel. The first 
assumption implies that the channel is modeled as ffequency-nonselective, while the 
second implies that the channel is slowly fading. The signal channel is modeled as a 
Ricean fading channel, hence, signal amplitude is a Riccan random variable [Refs. 2,8]. 
For Ricean fading, the total signal power consists of a direct signal component and a 
diffuse signal component, and the strength of the fading channel is characterized by the 
ratio of the direct signal component power to the diffuse signal component power. 
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The symbol rate is R^. For MFSK with M order modulation, the corresponding bit 



rate is = log 2 (M) . For L hops per symbol, the hop rate is = LRs- The spread 
spectrum bandwidth, W, is considered very large compared to the hop rate. 

B. PARTIAL-BAND NOISE JAMMING ANALYSIS 
1. Problem Development 

The partial-band noise interference when present is assumed to be in all branches 
of the MFSK demodulator and affects each chip of the dehopped signal with probability y 
where y is the fraction of the spread bandwidth being jammed. Hence, the fraction of the 

spread bandwidth not being jammed is 1 - y. If the average power spectral density of 

the interference is Ni/2 over the entire spread bandwidth, then the power spectral 
density of the partial-band interference when present is N^/ly. The power spectral density 

of the thermal noise and other wideband interferences which are modeled as additive 
white Gaussian noise is Ny2. Consequently, total noise power spectral density is N„/2 in 
the absence of partial-band interference; otherwise. 



is the total noise power spectral density when narrowband interference is present. 
If the equivalent noise bandwidth of each detector branch in the MFSK 



with probability 1-y when no jamming is present. When jamming is present, the total 




( 6 ) 



demodulator is B Hz, then the noise power received in a given hop is aj = =Noi5 



noise power in a given hop k is = + rs\ = 



0^1 
V y 



+ No Ji5 with a probability of y. 
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We assume that each receiver hop slot has the same noise equivalent bandwidth. The 



noise equivalent bandwidth of the receiver investigated in this thesis is B=Rf,. 

2 . Probability of Bit Error 

When partial-band interference is present the probability of symbol error for a 
MFSK receiver is 



where PJ^i) represents the conditional probability of a symbol error given that / of L hops 
of a symbol are jammed. Since each signal branch of the receiver is symmetric with the 
other branches, we can determine the probability of a symbol error , Pj , by considering 
the signal to be present only in branch one of the MFSK demodulator. The outputs of the 
other branches are assumed identical and independent (iid). 

For orthogonal MFSK the probability of a bit error is related to the symbol error 



The energy per bit as a function of the symbol energy and the modulation order is 




(7) 



by 




( 8 ) 



logjCAf) 



(9) 
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3. Probability of Symbol Error Under Partial-band Noise Interference 



Assuming the signal is present in branch one of the MFSK demodulator allows 
us to write the probability of a symbol error based on the conditional probability density 
functions that i of L hops are jammed.The conditional probability density function for the 

output of the branch containing the signalis fx^ixi I/) where X, is the random variable 
that represents the output of the signal branch. The conditional probability density 

functions for the non-signal branches aTQfx„(xm\i),m = 2,3,4...,M where the X„'s are 
the identically distributed random variables that represent the output of the branches that 

do not contain the signal. [Ref 3] The conditional probability of symbol error is 



^.(0 = 1 - 




( 10 ) 



for aWm^ 1 . 

Since the partial band interference may or may not be present in a hop, we must 
be able to differentiate between the two possibilities. Let the subscript n= 1, 2 denote that 
hop A: of a symbol has interference and has no interference, respectively. The diversity 
summer acts to add together all L independent hops in its branch. Of those hops / of them 
are jammed; hence, we can express the random veiriable at the output of the diversity 
summer (10) as 



Xn, = txmk„ = tx,„k,+ j: W=1.2,3...,M (11) 

/r=l k=\ k=L-l 
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a. Probability Density Function for the Decision Variable X„ 

The probability density function of the independent, identically distributed 
(iid) random variables X„,,, , m= 2,3,...,M that represent the demodulator branch outputs 
not containing the signal for hop ^ of a symbol is 

r \ 



f^rnki^mk) — 






-Xmk 

Vial ) 



u(xmk) m = 2,2,...,M 



( 12 ) 



where M(»)is the unit step function. [Ref. 9] 
The Laplace transform of (12) is 



Fxjs) = 



1 

2ct^5+ 1 



(13) 



Since each hop is independent, from (1 1), we can express the Laplace transform of the 
conditional density function,yA-„(X/„l/) , as 

^x„(^l0 = (^Io]‘ * (14) 



Direct inversion of (14) leads to an infinite series of confluent hypergeometric functions, 
but this proved difficult to program and slow to execute. To evaluate (34) we first invert 
the individual portions consisting of either all hops jammed or all hops not jammed. The 
conditional density function for the decision vtiriable is then obtained by convolution. 

The Laplace transform pairs are 




X'-' 

(/-I)! 



exp 




(15) 



and 
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(16) 









Ua^2>' 



L-i 



^L-i-\ 



(I-/-1)! 



exp 



r \ 

-X 

Ua^2>' 



w(x) 



where and < 3^2 represent the noise power in a hop experiencing interference and a hop 

experiencing only thermal noise, respectively. To recover fx„ixm) , the two m-Erlang 
random variables above are convolved . Hence, the inverse Laplace transform of (14) is 

obtained from a convolution of the individual right hand sides (RHS) of (1 5) and (1 6) as 



dt{\l) 





( 1 


( \ 
-t 


(x™-/)^-'-' 




UaJ 


[icl,] (/-!)! 


^2aii j 


(L-/-1)! 




^ 2a^2 7 



Applying the binomial theorem to the integrand of (17) and defining 



X = 






'^2a^i J 



V2a^2^ 



1 



(L-i- !)!(/-!)! 



(18) 



we can express (17) as 



l-l-x ( L-i- \ ^ 



^0 



xU-0"-'-'-*exp 



-t 



^ 2aJi 2a^2 J 



dt (19) 



Interchanging the order of integration and summation and integrating term by term, we 
get [Ref. 1 0, equation 3.351.1] 



/v„(xmli) = ^ £ 

h=0 \ 



where for brevity 



r T • 1 ^ 

L-i-\ 



. r-x. (-i)(^+/-i)!xi"'-'-' 

(- 1 ) ' *expi I Z 7177^20) 



a = 



2a 



( 21 ) 



k\ 



and 
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(22) 



P = 



1 



The special cases of all hops jammed {i=L) and no hops Jammed (/=0) are 

obtained directly from (15) and (16) respectively, since the analysis leading to (20) is not 

valid in theses cases. The remaining task is to evaluate the cumulative distribution 

function of y^„(x„|/) that is required for computing the probability of symbol error, (10). 

The computation required for the non-signal branch contribution to the probability of 
symbol error is fx„(xm)- There are two avenues to compute this, the first based on 



direct integration or, equivalently, the second based on an inverse Laplace transform of 
Fx„(s\i)ls. 



Directly integrating (20), valid when 0 < / < Z, we have 



f 

Jo 



Z.-I-1 



' r l-\j (Ct')P^ ' ^ ^ 

fx.(Frn\i)dXm- I ( D 



/z-/-n 



x(Z + /-l)! 



{L-i-\-k)\ 

(a-P)'"' 



1 

pi -,-1 



L-i-\-k L-\-\-k-p 

exp(-px^) I 

/7=o (L-i-l 



- ‘£‘ ,, ^ - exp (-ox.)T , 

1=0 (^ + /-l-/)(a-p) [cc 9=0 (Z-2-/-^)!a^+' J _ 

This equation is simplified by expanding the binomial coefficient and canceling like 
terms to obtain 



(23) 
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1 (A:+/-l)! f a'P* 

(/-I)! h k\ [(-l)'(p-a) 



l-exp(-px;„) S 

p=0 



P'- 



P^~' (Z,-2-/)!(-l)^~'~' 

(L-i-l-k)\ 1=0 (^+/_ 1 



1 -exp(-ca:;;,) H 

9=0 P'- 



(24) 



when 0 <i <L. For the case of all hops jammed, / = L, the cumulative density function is 









2a^ 



‘ -expi^iiic 



/H^ 






2a 



/H 



= 1 - exp • 






A* 



2a^i J *=o A:! 1,20^, J 



and for the case of no hops jammed, / = 0, it is 

fjfx„(x)dx= 



/ \ Z- 

1 1 u-1 



KlojiJ 



(i-1)! 



exp' 



2a 



>dx 



kl. 



(25) 



= 1 - exp • 



■sA 



/ \ 



2a*2j^o k\\2al2J 



A final useful form is obtained by making the linear transformation of the random 
variable X „, . Recognizing that a = — ^ and P = — ^,then 



2a 



2a 



(26) 



a _ 2aj^ _2l - 

P ■ -4 " a?" 

2^2 ' 

Then making the change of variables y = Pxm in (24) we get 



( 27 ) 
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(NtY 



O'-l)! k=o n\ [(_i)'(7Vt-1) 



k+i 



l-exp(->^) Z 



p=0 P 



(Z-0'- io (Z+/-i-/)!(iV7.-l)'+' 



, f 



(28) 



which is useful for numerical integration purposes. Applying the same substitution in (25) 
and (26), we get 



L-l , 

Fy(v\L) = 1 - exp {-y} S ^(v) 



(29) 



and 



Fr(ylO) = 1 - exp 



iLlyll'iLV 

Nrlto k\ ^Nt) 



(30) 



for the case of all hops jammed and no hops jammed, respectively. 

Previous efforts to determine the cumulative distribution for the non-signal 
branches, the second technique, is based on a numerical inverse Laplace transform of the 
function j x Fx„(s\i) [Refs. 3, 4]. Despite the finite sums of exponential-like terms in 
(28), this second method proved more accurate for large values of L . 

b. Probability Density Function for the Decision Variable X, 

Continuing with the assumption that the signal is present in branch one of the 
MFSK demodulator accompanied by Gaussian noise , the random variable 2T, is described 
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by the conditional density function/v, (xi I/) . Given a signal amplitude of Jl Ok and 



before diversity combining, the probability density function of X , is 



/x,,(xij«*) = -^exp 

2(3 k 






x\k + 2a\ 

2al 



lo 



cikj2x\k 






u(xik) 



(31) 



where /o(»)is the modified Bessel function of the first kind and order zero. The channel 
is assumed to have Ricean fading; therefore, a* is a Ricean random variable with a 

probability density function given by 



= -^ exp 



-{al+o} 

2a^ 



j ( Oka 



u{ak) 



(32) 



where is the average power of the direct signal component and 2g~ is the average 
power of the diffuse signal component [Refs. 2, 11]. The total average signal power of a 



hop kofa symbol is + 2 ct^ and is assumed to remain constant from hop to hop. 

The random variable X\k when hop ^ of a symbol has interference and no 
interference is denoted by X\k^ and Xuj respectively. The conditioning is removed by 

integrating the product of (31) and (32) with respect to a* from zero lo infinity with the 
result 



/Vi*„(^U") 



1 

2(^jL + 2cj2) 



exp 



-1 


( o'*' 

xiA„ + 2a- 




/(. 




2 


1 + 2a- j 




+ 2a- j 



(33) 
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for x\k„ ^ 0 where implies a jammed hop and implies no jamming. Since each 
hop is independent, from (1 1) we see that the Laplace transform of the conditional 

density function for the signal branch,yA',(xi I/) , is 

(34) 



The Laplace transform of the density function describing a single hop, fx^i^^xk) , is 






^ Pat„ 



(35) 



where 



pA:„ — 



>{2a^ + ol) 



(36) 



When(35) to the c„ power and the inverse Laplace transform is taken, the result is 
equivalent to taking the c„ fold convolution, denoted by . Letting c, = / and 
Cj = L - /, so that elements subscripted with a 1 correspond to jammed hops while 
elements subscripted with a 2 correspond to hops experiencing only thermal noise, then 






(2c„a^)‘ 



xexp[-p*„(xu„ +2c„a‘)] x /£.„.|(2p*„.y2c„aLtu, (37) 

where lc„-x represents the modified Bessel function of integer order c,-l . (Refs. 3,11] 
Unlike the case of the non-signal branch analysis, no analytic solution for the final 
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convolution for the decision variable Xi has been found. Symbolically, the probability 
density function for A", at the output of the diversity summer is described by 



fxi(xi lO = [/x„_ (XU, I/)]®' 0 



(38) 



C. NUMERICAL EVALUATION OF PARTIAL-BAND NOISE JAMMING 

For levels of diversity up to L- 20, the conditional probability of symbol error is first 
computed based on the signal energy-to-thermal noise power spectral density ratio, 
signal-to-jammer power spectral density ratio, signal fading direct-to-diffiise ratio, and 
fraction of spread spectrum bandwidth jammed. These input parameters form the basis 
for the probability of symbol error which is converted to probability of bit error using 
(8). To summarize 




cf-k _ direct signal power/hop 
noise power/hop 



(39) 




noise power/hop 



(40) 



Eb average energy per bit 



(41) 



No thermal noise power spectral density 



Eh average energy per hop 



(42) 



Nj interference power spectral density 




(43) 



and 0.0 < Y < 1.0 . 
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The computation of the probability of symbol error is achieved for each L by 
weighting and summing the probabilities of symbol error conditioned on i of L hops 
jammed. A numeric integration is performed over the product of the signal path 
contribution and the non-signal path contribution (10). 

1. Non-Signal Path Contribution 

The cumulative distribution function is computed using the analytic expression of 

finite summations given by (28) for 0 < / < T for small L. However, it is computationally 

faster and results in less round off error for large values of Z, if a numerical inverse 

Laplace transform is performed on j x F;k-„(5l/). The special cases of all hops jammed 
(/=Z) or no hops jammed (/=0) can be computed directly from (25) or (26), respectively. 

The convergence of the inverse Laplace transform is accelerated by the Euler 

transformation [Ref. 12] usually taking about 60 terms to reach relative errors on the 

order of 1 0"’. However, certain combinations of input parameters results in slower 

convergence. The inverse Laplace transform algorithm is limited to not more than 1000 

iterations. 

2. Signal Path Contribution 

For the special cases either of all hops jammed or no hops jammed, the required 

probability function fx^{x\ If) is given by (37). In the case of all hops jammed n=l, 
c, = L, and Cj =0, and in the case of no hops jammed n=2, c, = 0, and Cj =L. In both of 

these cases, (10) can be solved analytically, but the results are so complex that numerical 

evaluation is easier and more straightforward. When and i Z,(38) is evaluated by 
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the numerical inversion of (34) after which the probability of symbol error, (10), is 
evaluated numerically. 

The partial-band jamming fraction that yields the worst-case performance for the 

conventional receiver is obtained experimentally by computing the probability of bit error 

F F 

as a function of y for fixed values of and To cover the broad range of channel 

No Nj 

fading severity, results are obtained for several values of the direct-to-diffuse ratio. For 
weakly fading channels, exhibiting a strong direct signal component, a direct-to-diffuse 
ratio of 100 is used. This is essentially a Gaussian channel. For strongly fading channels, 
considered as nearly Rayleigh channels, a direct-to-diffuse ratio of one is used. Typical of 
a moderate Ricean fading channel, a direct-to-diffuse ratio of ten is used . 
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D. PARTIAL-BAND NOISE JAMMING RESULTS 



To determine the worst case partial-band jamming performance of the FFH/MFSK 
receiver the probability of bit error is computed as a function of y. The process is repeated 
for levels of diversity up to twenty, which represents the upper end for practical systems 
in today's technology. Figure 2 through Figure 9 display the probability of bit error for 
various levels of signal- to-noise ratio, total jammer power, and modulation order. The 
figures are calculated with the assumption that the signal direct-to-diffuse ratio (RJ is 
constant during a hop duration. Figure 2 demonstrates the performance for moderate 
signal-to-noise ratio with a near Gaussian channel. The signal-to-thermal noise ratio, 
Eh/Nj, , is the signal energy contained in one hop for a conventional BFSK receiver. The 
table below summarizes the cases displayed in Figure 2 through Figure 9. 



TABLE 1. PARAMETERS FOR WORST CASE ANALYSIS 



FIGURE 


Eh/N„dB 


EMdB 


Rfe / Ricean Fading 


Figure 2 


13.35 


3 


100 /Very Weak 


Figure 4 


13.35 


3 


10 / Moderate 


Figure 4 


13.35 


10 


100 /Very Weak 


Figure 5 


13.35 


10 


10/ Moderate 


Figure 4 


10 


3 


10 / Moderate 


Figure 7 


16 


10 


1 0 / Moderate 


Figure 8 


20 


3 


1 / Strong 


Figure 9 


10 


0 


1 / Strong 
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Pb 

0.300 

0.200 

0.100 



0.050 
0.030 

The surprising result is that despite increased levels of diversity, the performance 
does not greatly improve. There will be more to say about this later, but the same trend is 
observed in Figure 3 where there is Ricean fading of the signal . 




1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 

Number of Hops 



Figure 2. Worst Case Partial-band Jamming 
Eh/N,=13.35dB, EJNf=3dB, R^^lOO 
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Number of Hops 



Figures. Worst Case Partial-bandJamming 
Eh/N =13.35dB, Eh/A0=3dB, Rb=10 



Another, surprising result is that increasing the modulation order to four or eight 
does not provide much improvement. Figure 2 and Figure 3 display pessimistic 
performance from the communicator's viewpoint. However, they assume a high level of 
jammer power. We next investigate the performance when the jammer power is just one 
tenth of the communicator's power when Z,=l ( Figure 4 and Figure 5) . Since we are 
assuming 1000 hop slots, this is equivalent to assuming a signal-to-jammer power ratio of 
-20dB for BFSK with no frequency-hopping 
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Figure 4. Worst Case Partial-band Jamming 
Eh/N =13.35dB, Eh/AO=10dB, Rb=100 

Pb 




1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 

Number of Hops 



Figure 5. Worst Case Partial-band Jamming 
Eh/N<,=13.35dB, Eh/AO=10dB, Rb=10 
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A comparison of Figure 4 with Figure 5 illustrate the influences of Ricean fading 
on the probability of bit error when the communicator enjoys a significant power 
advantage over the jammer. Next we consider a situation where the communicator has a 
larger initial signal -to-thermal noise ratio, i.e. 16dB, to see if at higher levels the jammer 
could mitigate the effects of diversity with partial-band noise jamming. Figure 6 and 
Figure 7 display these results. 
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0.150 

0.100 
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Figure 6. Worst Case Partial-band Jamming 
Eh/No=16dB, Eh/A0=3dB, Rb=10 
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Pb 




Number of Hops 



Figure 7. Worst Case Partial-band Noise Jamming 
Eh/N„=16dB, Eh/AO=10dB, Rb=10 

To examine the extremes of channel fading on the received signal, the case of 
near Rayleigh fading is examined with a direct-to-diffuse ratio of one. Since fading is 
severe, we would expect poor performance even without jamming, but with some 
improvement added by the increased diversity. The signal energy per hop is examined at 
20 dB (Figure 8) and at lOdB (Figure 9). 
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Figure 8. Worst Case Partial-band Jamming 
Eh/N =20dB, Eh/A0=3dB, R,=l 



P b 




Number of Hops 



Figure 9. Worst Case Partial-band Jamming 
E,/N =10dB, E^/AO-OdB, R,=\ 
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In both Figure 8 and Figure 9 the improvement offered by diversity initially 
lowers the probability of bit error, but not to the extent that it can overcome the errors 
introduced by partial-band jamming. It should be noted also that the high 
signal-to-thermal noise ratio is received after passing through the near Rayleigh channel. 
Therefore, it has a diffuse signal component equal to the direct signal component. This 
represents a very poor channel for the communicator that is further degraded by 
partial-band noise jamming. 

The conclusion to be drawn from the preceding series of performance plots is that 
for a wide range of signal-to-thermal noise level, jammer power, and channel conditions, 
the increase in diversity in the FFH/MFSK receiver alone is not enough to overcome the 
effects of partial-band noise jamming acting in concert with multipath effects. The 
critical parameter then is the fraction of the spread spectrum bandwidth, y, over which the 
jammer power is spread. The Figures 10 through 17 demonstrate the fraction of jammed 
bandwidth that provides the worst case probability of bit error. 
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Figure 10. Fraction of Jammed Bandwidth 
Eh/N„=13.35dB, Eh/iV;=3dB, Rb=100 
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Figure 11. Fraction of Bandwidth Jammed 
Eh/No=13.35dB, Eh/A0=3dB, Rb=10 
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Figure 12. Fraction of Bandwidth Jammed 
Eh/N =13.35dB, Eh/A^^lOdB, Rb=100 
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Figure 13. Fraction of Bandwidth Jammed 
Eh/N„=13.35dB, EJNj=\0dB, Rb=10 
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Figure 14. Fraction of Bandwidth Jammed 
Eh/N =16dB, Eh/A^^3dB, Rb=10 
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Figure 15. Fraction of Bandwidth Jammed 
Eh/N„=16dB, Eh/A0=3dB, Rb=10 
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Figure 16. Fraction of Bandwidth Jammed 
Eh/N,=20dB, Eh/A>=3dB, Rb=l 
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Figure 17. Fraction of Bandwidth Jammed 
EH/N=10dB, Eh/A0=OdB, R,=l 
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The worst case values y are determined by computing the probability of bit error 
as a function of y. Figure 10 through Figure 1 7 show that , in general y, depends on the 
hop energy-to-thermal noise ratio (Eh/No ), the hop energy-to-jammer power spectral 
density (Eh/A^), the level of diversity {L\ the modulation order {M), and the channel 
fading direct-to-diffuse ratio. Further, the shape of the worst case y is generally the 
same, showing an initial drop as L increases with a gradually lowering slope. It is 
important to note that the lowest levels of y are still above 10'^ which represent only one 
jammed hop slot when there are 1000 hops in the frequency-hopping system. 

That the increase in diversity is unable to mitigate the partial band jamming is 
also indicated by looking at the trend of performance for fixed values of y. Figure 10 
through Figure 17 show that performance improves linearly as diversity increases. This 
is not unexpected since the constant energy per hop system increases the total energy per 
bit as diversity increases, however, the change in slope of the performance curves is 
indicative of the effect contributing to the worst case analysis. Figure 18 through 
Figure 23 display the performance for fixed values of y. 
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Figure 18. Performance for fixed values of gamma 
Eh/N„=13.35dB, Eh/AO=10dB, Rb=10 
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Figure 18 displays the performance in a moderate Ricean fading channel where 
communicator enjoys a fair signal-to-jammer power advantage. The performance 
improves linearly as diversity increases. However, the slope of the performance curve is 
significantly reduced for decreasing values of y. Further, we can see that for a given L 
there corresponds a worst case value for y which becomes less sensitive as diversity 
increases. Figure 19 displays the performcuice when the signal enjoys less power 
advantage relative to the jammer. The trend is that for lower signal-to-jammer power 
ratio performance improves more slowly. 
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Figure 19. Performance for fixed values of gamma 
Eh/N„=13.35dB, Eh/A0=3dB, Rb=10 



The worst case fraction of bandwidth jammed is influenced by channel fading 
(Figure 10 through Figure 17). Likewise, the slope of the performance curves is 
influenced by channel fading . As the fading becomes more severe, i.e., near Rayleigh, 
the performance improves more slowly as demonstrated in Figure 21 through Figure 23. 
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Figure 20. Performance for fixed values of gamma 
Eh/N„=13.35dB, Eh/A^^lOdB, Rb=l 
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Figure 21 . Performance for fixed values of gamma 
Eh/N„=13.35dB, Eh/A0=3dB, R,=l 
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Figure 22. Performance for fixed values of gamma 
Eh/N =13.35dB, Ey,/Nf=\0dB, Rb=100 
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Figure 23. Performance for fixed values of gamma 
Eh/N„=13.35dB, Eh/A0=3dB, Rb=100 
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III. MULTI-TONE INTERFERENCE 



A. SYSTEM DESCRIPTION 

In addition to partial-band noise jamming, spread spectrum communications can 
experience tonal or narrowband interference in one or more hop slots. We can similarly 
expect that as spread spectrum bandwidth grows the likelihood of tonal interference will 
also grow. This interference may arise from hostile sources as well as non-hostile ground 
communications, satellite transponders, and radars. The multi-tone interference problem 
is similar to the partial-band noise problem. 

This chapter analyzes the probability of bit error arising from an "intelligent Jammer" 
who knows the spread spectrum bandwidth, fV, and the modulation order. This 
intelligent jammer is able to place a single interfering tone in one or more hop slots, but 
has no knowledge of the pseudonoise (PN) sequence driving the hopping pattern. This 
represents a very intelligent jammer. We then relax this situation and compare it with the 
performance of the jammer who places an interfering tone in none, one, or both of the 
symbol tone locations in a particular hop slot. 

It is likely that the interfering tones may experience multi-path effects in a manner 
similar to the signal tone. However, since the interfering tones and signal tone can arrive 
from different paths, the interfering tones will in general arrive at the receiver with a 
different direct-to-diffuse power ratio from that of the signal. In this thesis, the 
performance when both the signal and the jammer experience fading is evaluated. 
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B. MULTI-TONE INTERFERENCE ANALYSIS 



1. Problem Development 

First we consider the case where the jammer places at most one tone in each fast 
frequency-hop slot. Only the worst case situation where the jamming tones exactly 
coincide with frequency-hopped symbol tones is considered. The situation may be 
visualized with the aid Figure 24. 

FFH/BFSK 
N slots 




W- 




|C| tones jammed 



Figure 24. Multi-tone interference in FFH/BFSK 
The FFH/BFSK receiver is a special case of the FFH/MFSK receiver described 

earlier. The orthogonal design is selected such that is the system bit rate and each 

signal tone is placed b apart. In Figure 24 the bandwidth of a single hop. B, is 2.^6 b • 

The spread spectrum encompasses W Hz and is divided into N hop slots. Although 

numbered consecutively for illustration, the pseudonoise modulation will create a 
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random-like sequence to the hopping pattern. In this thesis, iV=1000 is choosen as a 
representative number of hop slots. Further, the number of interfering tones, q, must be 
an integer 

\<q<N (44) 



Since the jammer spreads his available power, J, equally among q tones, each 
interfering tone has a power of J^=J/q. Earlier work focuses on the simplifying case 
of no thermal noise [Ref. 13]. In the no thermal noise case an error can only occur if the 
hop is jammed and the the interfering tone power is larger than the signal power. In this 
instance the jammer can make most efficient use of his power by choosing the number of 
tones to jam as 



q = INTO 



J_ 

Pc 



(45) 



where is the carrier power of the desired symbol and INTG[ ] represents the greatest 
integer less than the argument. [Ref 13, pp.596-598] In this thesis, the probability of bit 
error is obtained using fixed values of q to determine a worst case scenario which 
includes the effects of both signal and jammer fading as well as thermal noise. The 
approach to this problem is to first detemine the probability of bit error for a conventional 
BFSK system without frequency-hopping and then generalize those results to the 
FFH/BFSK system. 
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2. BFSK Analysis with Single Tone Interference 



For the binary case we begin with the familiar expression for the probability of 
bit error for the noncoherent FSK receiver 

Pb=Psx >/?2l52] (46) 

where i?, and R 2 represent the demodulator outputs of the branch detecting signal 5 , and 
S 2 , respectively. Making the assumption of equally likely signaling, we have 

=5 (47) 

and 

Pb = \[PMl >R^)^Pr{R^> /? 2 )] (48) 

The received energy consists of signal energy, thermal noise, and, if jammed, the 
interference signal 

r{t) = Si{t) + n{t) + sj{f), (49) 

where is the ith signal tone, Sj is the interference tone, and n(t) is additive white 
Gaussian noise with power spectral density N,, / 2. Since the receiver structure is 
symmetric, we can simplify the analysis problem by considering the case where branch 
two is jammed. 

First we examine the case where the signal and interference occupy different 
frequencies. Then 

Pr(error\si) = Pr(R2 > 1^1) (50) 
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We express the received waveform as 



r(0=^/2fl cCOs((Oi/ + 0) + Jl ajcos(( 02 t + Qj) + n(f) 



(51) 



where and Oj are modeled as Ricean random variables. As previously discussed the 
received amplitudes and phases of each tone is a random variable due to multipath 
effects. For the noncoherent detector design, the output of each branch can be expressed 
as a non-central Chi square random variable [Ref. 2]. For the signal branch 






2cj" 






2a' 







do 


L C7^ J 



u(ri) 



(52) 



and for branch two containing the interference tone 



fR,(r2\aj,si) = ^exp 
IG 



{r2+2alTl) 


> Ir. 




2g^ 




2 V 

i G^ J 



u(r2) (53) 



where a = The probability of error (50) is 



PriRi > R\\ac,ajS\)= ^fR,R2(r\.r2\aj.ac.s\)dr2dri (54) 



This represents the probability conditioned on the ^lmplitudcs of a. and Oj . The 
orthogonal receiver structure and independent fading of the interference tone and the 
signal tone allow the joint probability density function to be separated. Hence. (54) 
simplifies to 

PriRi >R\\si) = J^/«,(rilflc,-yi)[j^^//?2(A-2l<37,-!i'i)a'/'2]i/''i (55) 
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The evaluation of (55) is simplified by applying the quadratic transformation to convert 

(52) and (53) to Ricean random variables [Ref. 11]. If ri = , then = 2vi . Using 

dv\ 

this substitution in (52), we obtain 



/r,(vilafc,^i) = ::p7exp 
2a^ 



(v^ +2alfl) 



2ct^ 



Ko 



JlacTbV] 



m(vi) (56) 



Similarly for the branch output containing the interference tone 



fv^(v2\aj,Si)=^exp' - 

2<t^ 






2a^ 



h 



JlajTbVi 



w(v2) (57) 



Now (55) can be expressed as 



Pr(V2 > F|l5i)= J^/) ,(V|bc,5l)[j^ /l':(V2la7.^l)<^^2 



dv\ 



(58) 



The inner integral equation has no simple solution, but is represented symbolically as 
Marcum's Q function which is defined as [Ref 1 4] 



QiP-^ P) = C vexp + o}) /o[av]c/i 



(59) 



Hence, (58) is expressed as 



/’r(v 2 > v,bi)= Q^JlajTblo,'^') x ^expj- ^' 



JlucTb 



V| 



dv\ (60) 



which can be evaluated to obtain [Ref 15] 



> V, I. ,) 4[ 1 - e(2#. 2#) . 2#) 



( 61 ) 
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Equation (58) represents the contribution to the probability of bit error when the signal 
and jammer occupy different tones within a hop. It remains to determine the contribution 
due to a collocation of interfering and signal tones. That is 

Pr(error I ^■ 2 ) = Pr{R 1 > I 52 ) (62) 



where 



A(r,) = ^expj^ 



u(r\) 



(63) 



Since the signal and the interference are in the same tone location, the output of 
the in-phase and quadrature portions of branch two due to signal and jamming tones is 



JlTb [flccos (0) + ajcos ( 87 ) (64) 

and 

^2 Ti [ricsin (0) + ojsin ( 87 )] (65) 



respectively. The phase of the jammer is generally different from that of the signal. 
Squaring and summing to form the branch output produces 

i?2(0 - 87) = 2T^[ac + ^*7 + 2flc07(cos 0 cos 07 + sin 0 sin 87] (66) 

which is simplified to 

R2{h) = 2T^^al+a] + 2acajcosb\ (67) 



where5=0-0j. Now, replacing 2a/r,,‘ with (67) in (53) we get 



/« 2(''2 Iflc, «7, 5, 52 ) = exp 

2a^ 



{t 2 + 2Tl{al j 



2a' 



No 



^j2r2 (0^ + 07) 



u{r2) ( 68 ) 



Now, 
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Pr(Ri >R2 \s 2,?>)= ^^fR 2 {r 2 \ac,aj,h,S 2 )y-^^^jR,(r\\s 2 )dr\'^dr 2 



(69) 



Thus, the probability of error arising from collocation of the jammer and the signal is 
computed from a probability density function conditional on 5. In this case the 
non-signal branch is unaffected by the jamming tone. Substituting (63) and (68) into 
(69), we get 



Pr{R\ >i?2l5'2,5) = ^xexp 



+ lacQjCOS 5 ^ 



(70) 



This can be simplified into familiar terms by replacing 



in ^ 2Tt ^Et ^ 



(71) 



and similarly for the jammer power term 



„i ^ b 
Oj ^ 






in 

No 



No 



yj 



(72) 



Then the conditional expression in (70) then becomes 



Pr(R 1 > /? 2 1^ 1 , 6) = (i) exp j -^(y* + yj + 2 Jjtjj cos 6 



(73) 



The nuisance parameter 5 may be removed by multiplying by its probability density 
function and integrating the conditional density. That is 



Pr(R\ >R 2 \s\) = -^ f -y{yb+yj + 2Jybyj cos 5 



d8 



(74) 
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where 6 is modeled as a uniform random variable. Since the integral is only over 5, then 
[Ref. 10] 



PriR\ >./? 2 li’i) = exp jY(yi + y7> exp { ^y^ cosd} (75) 



which can be evaluated to obtain 



Pr{R\ >./?2l^i)= (^) xexp' ~(yfc+yj)J x/o[yy^] (76) 



Therefore, the phase relationship between the signal and the jammer plays a significant 
role as evidenced by the modified Bessel function term. This makes physical sense. 
When the phase is the same or nearly so there is reinforcement of the signal tone by the 
interfering tone, but when the phase difference approaches 1 80 degrees there is 
destructive interference. 

This allows us to compute the probability of bit error in the absence of fading by 
substituting (61) and (76) into (48). 



Pb(oc,aj) 




+exp- 



zl± 

2a^ 



(«c + I -^0 



acOjfl 



(77) 



An equivalent expression using (71) and (72) is 



7’/.(y*,yy) = 



+Q{jyJ,Jyb) + exp 



y(Y6+Yy) 






( 78 ) 



Using a Q function property [Ref 15, pp. 396], we can reduce this to an expression 
containing just one Q function as 
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(79) 



Pb{ac,aj) = \-\ X Q 



ajTb OcTb 



lo 



acajTy. 



-1 



X exp • 



zli 

2a^ 



cic+aj^ 



Now, the amplitudes of the received interference tone and signal tone are Ricean random 
variables where 
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Where in each case the second moment can be expressed as 

9 9 9 

a, =aj +2a, fori=J,c 
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u(aj) 



u(ac) 



(80) 



(81) 



(82) 



We first remove the conditioning on in (79). The technique is to take each piece of 
(79), multiply by (81), and integrate over all possible values of a^. Hence, for the first 
piece 





fAciflc)dcic ~ 
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^iTl 
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The integral of the product of the second piece and (81) is 



( 84 ) 
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Making the substitutions 
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( 86 ) 
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in (85), we can use [Ref. 10, equation 6.633.4] to obtain 
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for the second piece of (79). The remaining piece to be determined is 

ClcCLc 






f 2 . 2 

' + a^ 



2 a^ 



■/o 



da 



(87) 



( 88 ) 



(89) 
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Using [Ref. 10, equation 6.631.4], we evaluate (89) to obtain 
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(90) 



Finally, combining (84), (88), and (90) with (79), we get the probability of bit error when 
the signal experiences Ricean fading as 
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(91) 



which is now conditional on Oj. Expressing (86) in terms of 



C/.= 
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No 
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( 93 ) 



and the direct-to-diffuse power ratio 
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( 94 ) 



R =^ = 
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we get 
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Using [Ref. 15 ] 
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in ( 95 ), we obtain 
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( 97 ) 
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We now employ a similar approach to remove the conditioning on Qj. Integrating 
the product of (80) and (97) over all possible receiver jammer amplitudes, Oj, we get 
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(98) 



where 
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An alternate expression for the probability of bit error in the case of both signal and 
jammer fading in terms of weighted Q functions is 
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( 102 ) 



However, since there is no closed form solution for the Q function, in practice it is easier 
to compute the performance using (98). 
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a. Special Cases of BFSK with Single Tone Interference 

The final aspect of the analysis is to consider how the derived result 
simplifies for special cases, both as a check on earlier work and to see limiting cases. As 
check on (98), when there is no signal the probability of bit error reduces to one half as 
expected. 

Next consider the case of no jamming. In this case, §j= = 0, and the 

probability of bit error reduces to 



which is the usual BFSK result. 

Another limiting case is the performance when both the signal and the 
jammer experience Rayleigh fading. In this case both l^^, and are zero and all the 
received energy is in the respective diffuse component. In this case (98) reduces to 



This displays an inverse linear relationship similar to Rayleigh fading without jamming 
[Ref 2]. Performance is improved if we change the case slightly to that where the 
jammer suffers Rayleigh fading while the signal experiences Riccan fading. Then 



Here any direct signal component that reaches the receiver will serve to drive the 
probability of bit error lower in an exponential fashion. 




(103) 




(104) 




(105) 
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Perhaps the most optimistic performance from the communicator's point of 



view is when the signal has no fading, but the jammer suffers Rayleigh fading. In this 
case, ^ 0, and (98) reduces to 



On the other hand , the most pessimistic performance is likely to occur when 
the signal suffers Rayleigh fading, but the jammer has no fading. In this case, 

= 0, and (98) reduces to 




1 






exp 




(107) 



This performance will be very poor unless there is a relatively large amount of diffuse 
signal energy received. 

3. FFH/BFSK Analysis with Multi-Tone Interference 

We now turn our attention to the application of (98) to spread spectrum and the 
partial band multi-tone jamming scenario. Suppose N is the number of hop slots in the 
FFH/BFSK system and the jammer transmits q interfering tones. When we credit the 
jammer with only placing at most one tone per hop slot then 

r = (^) (108) 

is the fraction of spread spectrum bandwidth jammed where y is restricted to 

7-<y<1.0 (109) 

N 

Now, 
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Pb = I hop jammed] + {^^-j^jPr[error\hoY> not jammed] 



where 



Pr[^^^orlhop not jammed] = Pr\error\BFSK\ 





(110) 



To apply our earlier result for the BFSK receiver we must define our jammer parameters 
on a per hop basis. If is the average total jammer power, then the average total 
jammer power per hop is 



y Y Y • 



( 111 ) 



which implies and are replaced by N^j/q and Nl^j/q , respectively, in (98). Now (1 10) 
is given as 







( 112 ) 



In the case of FFH/BFSK with one hop per bit, no fading, and no thermal noise 
(N„=0), an analytic expression for the worst case number of jamming tones, , is 
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1 <INTGf^'i =qo<N 

\yb J 



(113) 



One aim of this thesis is to determine if this limited analytic result for remains true for 
FFH/BFSK over Ricean fading channels in the presence of thermal noise. This is 
accomplished by comparing the performance based on q^ and some offset values of q. 

4. FFH/BFSK Analysis Allowing Two Interference Tones 

We now extend the previous results to a situation which includes the possibility 
that both receiver branches contain an interfering tone. This is a relaxation of the 
intelligent jammer scenario employed thus far. In this situation, the probability of bit 
error may be expressed as 

Pb = Pr[no hops jammed] x P;.[errorlno hops jammed] 

+Pr[one tone jammed] x P,.[errorlone tone jammed] 

+/’;.[two tones jammed] x P,.[errorltwo tones jammed] (114) 

The first two conditional probabilities were previously determined in Section 2. 

However, the probability of those events occurring is a sample without replacement 
situation. That is 

Pr[one tone jammed] = P^[one tone jammed n one tone not jammed] 

Pr[0 tone jammed 1 1 tone not jammed] x Pr[\ tone jammed lO tone not jammed] 
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(115) 



V2A^A 2A^-ly' V2A^A2iV-lJ 

Similarly 

P.[no tones jammed] = (l “ ^) (l “ 
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I 2A^ A 2N-1 ) 



And lastly, the probability that both tones are jammed is 

Pr[ both jammed] = Pr[\ jammed n 0 jammed] 

_( ^ y 

V2A^A2A^-ly' 



(116) 



(117) 



All that remains is to determine the probability of error given that both branches are 
jammed. 

a. Extension of Single Interfering Tone Results 

Now that both branches of the BFSK receiver contain jamming tones, 
symmetry is restored . [Ref 11] Hence, we can determine the probability of bit error 
supposing that signal 5, is sent. An error occurs when the output of branch two is greater 
than the output of branch one and (54) yields 



/’r[eiTorltwo tones jammed] = 
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(118) 
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Alternatively, this may be expressed as 



Pr [error I two tones jammed] = Q 



ajTb 7 ’* Jqc +aj + lacojcos 6 
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(119) 



The conditioning on the received amplitudes of the signal and the jammer is removed by 
multiplying (119) with (80) and (81) and integrating over all values of and . In 
addition, the conditioning on the signal-to-jammer phase difference, 5, must also be 
removed. 

In order to make (1 19) valid for FFH, we substitute (111) into (119). Now, 
symbolically the probability of error in the case of both tones jammed, » is 
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( 120 ) 



61 



Substituting our previous results for BFSK and (120) into (114), we get the 



total probability of bit error 
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( 121 ) 



Further simplification results by using the substitutions 
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and 
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which implies 
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Then using (122) through (125) in (80) and (81) 



(125) 
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(126) 



and 
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exp 
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(127) 



b. Special Cases when Allowing Two Interference Tones 

As a consequence of the triple nested numerical integration required by (120) 
in order to evaluate (121), obtaining numerical results when two tones can jam a single 
hop slot is computationally intensive. Based on the results when thermal noise is 
neglected, we expect the performance will be better than when only a single interfering 
tone is allowed. 

When there is no signal fading 



Y7 = Ya and^A =0 (128) 

and when there is no jammer fading 

Yy = Yy = Oand^j = 0 (129) 

In the case of Rayleigh fading of the signal 

Ya = and (i;* = 0 (130) 

and 

Y7 = ^7andO = 0 (131) 



when the jammer is Rayleigh faded. Two limiting cases allow evaluation of the most 
pessimistic and most optimistic performance. Optimistically, from the communicator's 
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point of view, when there is no signal fading and Rayleigh fading of the jammer, the 
probability of bit error reduces to 



* y 2N 2N-\ 
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which requires two numerical integrations rather than three. 

The most pessimistic performance result from Rayleigh fading of the signal 
and no fading of the jammer. In this case. 
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Further simplification when neither the signal nor the jammer fade. In this 



instance (132) reduces to 
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(134) 



where a single finite numerical integration over the phase angle is required. In this thesis 
only (134) is evaluated. Further work in this area could consider the effects of Ricean 
fading in the both tones jammed case based on (120). 
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C. MULTI-TONE INTERFERENCE NUMERICAL PROCEDURE 



The crux of the numerical computation for this problem is the accurate and efficient 
computation of Marcum's Q function. Originally the Q function described a shorthand 
notation for the probability integral representing the output of a correlation detector 
containing a radar target in the presence of narrowband Gaussian noise [Ref 14]. 

The Q function, described by (59), has no closed form solution expressible in a finite 
number of terms. Frequently cited equivalent expressions containing an infinite series of 
Bessel functions, while valuable as analysis tools, do not ease the computational burden. 
The strategy is to test the input arguments for their magnitude, difference, and the 
presence of zeros. Based on the input arguments the program computes the value of the 
Q function in an appropriate subroutine and returns to the calling program. 

First, the presence of a zero passed as an input argument is tested. If present the 
following simplification results 



The second consideration is the magnitude and difference of the input arguments. If the 
product of the arguments is greater than ten and their difference is greater than five, then 
the Q function is computed using an asymptotic polynomial. If the magnitude of the 
product is greater than 1000, then the Q function is computed by directly integrating (59) 
using a large argument approximation for the Bessel function. For the semi-infinite 



0((x,O)=l.O 
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and 




(136) 
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integration required here, the Romberg technique is very efficient. The Romberg 
technique is a recursive implementation of Richardson's extrapolation which is continued 
until a user specified tolerance is reached. In theory and in practice it offers much faster 
convergence and smaller errors than the standard trapezoidal or Simpson technique. 

[Ref. 16] 

If none of these conditions is satisfied the program defaults to a computation based 
on a numerical inverse Laplace transform. The integrand of (59) is expressed in the 
Laplace domain as [Ref. 10, equation 6.643.4] 
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(137) 



Then the integration may be determined from the inverse transform of 



0(a,P)c^F(.s) = |x xexpl-(^)(l- 
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(138) 



This is useful for many intermediate values of the input arguments. 
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D. MULTI -TONE INTERFERENCE RESULTS 



First, the results for the scenario in which at most one tone is jammed per hop slot 
are presented. The performance for a wide variety of channel conditions is computed for 
several representative signal to thermal noise ratios. All figures are for one hop per bit 
(Z=l). For each channel case there are two prime questions to be answered. First, "How 
does performance vary for fixed values of qT Second, "Is the no fading analytic and no 
thermal noise worst case q still valid for fading channels?" To answer the second 
question the proposed worst case value of q, q„ ,is used to compute performance as well 
as some values of q offset from q^ . To show the trend, the results of q„+ 4 and q^ -4 are 
compared with q^ based performance. Table 2 provides an overview of the channel 
conditions considered. 
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TABLE 2. MULTI-TONE CHANNEL CONDITIONS 



NAME 


Eb/N„dB 


Rb/ 

Ricean Fading 


Rj / ! 

Ricean Fading | 


Figure 25 


13.35 


100 /Very Weak 


100 /Very Weak 


1 Figure 26 


13.35 


100 / Very Weak 


0.0 / Rayleigh 


Figure 27 


13.35 


0.0 / Rayleigh 


100 /Very Weak 


Figure 29 


13.35 


10 /Moderate 


10 / Moderate 


Figure 29 


13.35 


0.0/ Rayleigh 


0.0/ Rayleigh 


Figure 30 


13.35 


10 /Moderate 


1 / Strong 


Figure 3 1 


13.35 


50 / Weak 


50 / Weak 


Figure 32 


13.35 


25 / Low 


25 / Low 


Figure 33 


20 


10 /Moderate 


10/ Moderate 


' Figure 34 


20 


1 / Strong 


10 / Moderate 


Figure 35 


20 

1 


1 0 / Moderate 
1 1 


1 / Strong 1 



Initially, we consider the case where fading has very little influence to compare with 
a no fading analytic result. Figure 25 is an illustration of the performance obtained with 
a moderate signal-to-thermal noise ratio that typically provides a 10'^ bit error ratio in the 
absence of fading and interference. Thus, Figure 25 provides a good basis for 
comparison with the no fading case. The worst case performance is seen as the envelope 
of the fixed q curves. As the signal-to-jammer power increases, the performance quickly 
approaches 2x10 ^ which is very close to the no jamming performance. We can also see 
that the degradation inflicted by multi-tone jamming is strongly influenced the choice of a 
fixed number of interference tones. 
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Pb 




Figure 25. Eb/N„=13.35dB, Rb=100, Rj=100 
To bracket the range of expected performance the case of extreme of Rayleigh 

jammer fading and very weak signal fading is depicted in Figure 26. Here the 

performance is slightly improved at lower signal-to-jammer power ratios than is 

observed in Figure 25. Equivalently, the worst case envelope is slightly inclined more 

toward the y-axis indicating better performance. However, at low signal-to-jammer 

power levels, the performance is still relatively poor. 
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Eb/NJ (dB) 



Figure 26. Eb/N„=13.35dB, Rb=100, Rj=0.0 



71 



Pb 




Figure 27. Eb/No=13.35dB, R^,= 0, Rj=100 
At the other extreme, a pessimistic performance is obtained when the signal suffers 

Rayleigh fading, but the jammer has very weak fading. The jammer is able to reduce 

performance overall and maintain a high bit error ratio (BER) for higher signal-to-jammer 

power ratios. Even at Eb/AO=30dB, the effects of fading and jamming produce a BER 

near 4x10'^. 
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Figure 28. Eb/N„=13.35dB, Rb=10, Rj=10 
In Figure 28 we display the performance when both the signal and the jammer 

experience moderate Ricean fading. As expected the performance lies between the 

optimistic case of Figure 26 and the pessimistic case of Figure 27. In this case each tone 

contains the same direct-to-diffuse power ratio, but as Ei/Nj grows above OdB, the direct 

signal power is greater than the direct jammer power and is able to mitigate the influence 

of multi-tone jamming. 
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Figure 29. Ei/N„=13.35dB, Rb=0, Rj=0 

In Figure 29, the signal and jammer again share the same direct-to-diffuse power 
ratio, but now Rayleigh fading is assumed. The performance is worse overall compared 
with the moderate Ricean case and approximately 3dB worse at higher signal-to-jammer 
power levels. Also note that the worst case occurs at ^=1000. Therefore, the greatest 
degradation occurs when the jammer is forced into the multi-tone equivalent of 
broadband jamming. 
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Figure 30. Eb/N„=13.35dB, Rb=10, Rj=l 

In Figure 30, the channel conditions are more favorable to the communicator and 
represent a ten fold increase in the Rb over the R, from the performance depicted in 
Figure 28. At the critical value of OdB for the signal-to-jammer power ratio (EJNj) the 
performance in the favorable channel ( Figure 30) is 3.3x10"* while in the moderate 
Ricean channel (Figure 28) the same EJNj produces a BER of 2.1x1 0 \ So the fact that 
the jammer suffers more fading over the Ricean channel produces a small improvement 
for the communicator. 
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Figure 31. Eb/N„=13.35dB, Rb=50, Rj=50 
Figure 3 1 is similar to Figure 25, but includes less fading. The performance is 

essentially the same as the very weak fading performance, but with some slight 

improvement for ^Nj greater than OdB. As in earlier fading channel cases, the 

performance varies greatly for fixed values of q when the jammer has a power advantage. 



76 



Pb 




Figure 32. Eb/N„=13.35dB, Rb=25, Rj-25 
Figure 32 depicts another drop in the channel direct-to-diffuse ratio shared by the 

signal and jammer. At this level the difference in performance below the very weak case 

where the direct-to-diffuse ratio was 100 is noticeable. When the jammer enjoys a power 

advantage, the BER is still unacceptably high. The worst case envelope curve is nearly 

linear above Ei/No=0dB and nearly flat for all choices of q above OdB. 
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Figure 33. Eb/N =20dB, Rb=l 0, R,=l 0 

In Figure 33 we begin to look at the performance for larger signal-to-thermal noise 
ratios. An increase in Eb/No from 13.35dB (Figure 28) to 20dB (Figure 33) both in 
moderate Ricean fading, provides an overall improvement in performance. However, the 
improvement is significant when E^/NJ is greater than -12dB and provides nearly a factor 
of ten at OdB. 
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Figure 34. Eb/N„=20dB, Rb=10, Rj=10 

Figure 34 shows how rapidly performance can degrade even for a strong 20dB 
signal-to-thermal noise power ratio when the signal fading progresses from moderate 
(Figure 33) to strong Ricean fading. Similar to the situation observed in Figure 27 , the 
worst case performance is obtained when the jammer places an interfering tone in every 
hop slot. The jammer therefore, need not be very sophisticated to prevent efficient 
communications. In Figure 35 the roles are reversed. With the same signal-to-thermal 
noise ratio, the signal experiences moderate Ricean fading while the jammer experiences 
strong Ricean fading. 
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Figure 35. Eb/N„=20dB, Rb=10, Rj=l 

With moderate Ricean signal fading a particular value of is required to cause the 
worst case performance demanding greater jammer sophistication. The strongly faded 
jamming tone is much less effective when E^JNj is greater than OdB. Previously we 
observed an approximate tenfold drop in BER experienced for the E[/No=13.35dB case as 
the jammer fading grew from moderate (Rj=10) to strong (Rj=l) at Et/A(,=0dB. This 
trend is more pronounced at the Ei/No=20dB level as we transition from Figure 34 to 
Figure 35 where the improvement is greater than a factor of 25. 
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The influence of fading on the selection of q to cause worst case performance in the 
multi-tone environment remains to be answered. For FFH/BFSK in the absence of fading 
and thermal noise we expect the worst case to occur when q is chosen as the integer 
portion of the jammer- to-signal power ratio. A comparison of the following worst case 
performance curves will show them to be the envelope of the previous fixed q computed 
curves. 

Pb 




Figure 36. Worst Case Multi-Tone Jamming for 
Eb/N„ = 13.35dB, Rb=0, Rj=100 

The most pessimistic case for the communicator is pictured in Figure 36 where the 
signal has Rayleigh fading and the jammer enjoys essentially no fading. We see that in 
this situation the jammer does not need to be very sophisticated to hinder communication. 
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The results displayed in Figure 36 and in Figure 37 indicate there is essentially no 
difference in the performance and the +4 and q^ -4 cases. 

Pb 




Figure 37. Worst Case Multi-Tone Jamming 
Eb/N„=13.35dB, Rb=0, Rj=0 

The conclusion to draw from Figure 37 is the degradation in channel conditions for 
the jammer does not significantly assist the communicator when both are acting through 
Rayleigh channels. 
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For convenience we will denote the and q^-A collectively as Isq^. It is not until 
there is some appreciable direct signal power reaching the receiver that the q^ and 
begin to show some deviation. 

Pb 




Figure 38. Worst Case Multi-Tone Jamming 
Eb/No=20dB, Rfc=10, Rj=10, 

In Figure 38 the curves dip below the anticipated worst case curve. The trend 
progresses as the direct-to-diffuse ratio increases in Figure 39 and Figure 40 respectively. 
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Figure 39. Worst Case Multi-Tone Jamming 
Eb/N„=13.35dB, Rt=25, Rj=25. 

Since the value of q is required to be a positive integer greater than or equal to one, 
there is a small discontinuity in the q^-A curve. At those signal-to-jammer power values 
where the q^ value would be less than 5, the q^-A defaults to q=\ for the computation. 
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Figure 40. Worst Case Multi-Tone Jamming 
Eb/N„=13.35dB, Rb=50, R,=50, 



The salient observation in Figure 38, Figure 39, and Figure 40 is that the anticipated 
worst case performance in the no fading and no thermal noise case remains the worst case 
performance for the weakly and moderately faded channels with moderate thermal noise. 
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For comparison Figure 41, Figure 42, Figure 43 displays the composite worst case 



performance for various signal-to-thermal noise ratios. 

Pb 




Eb/NJ (dB) 



Figure 41 . Worst Case Multi-Tone Jamming for Several 
Values of Efc/No, Rt=l, Rj=l 

Figure 41 displays the performance when both the signal and the jammer suffer 
strong fading. The diffuse nature of the channel keeps the probability of bit error 
relatively high despite high signal -to-Jammer power ratios and serves more to the 
detriment to the communicator than the jammer. 
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When the channel has moderate vice strong Ricean fading the performance is that 
shown in Figure 42. Above -6dB EJNj there is a significant improvement in 
performance over that in Figure 41 . The noticeable change in slope for the 20dB curve is 
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Figure 42. Worst Case Multi-Tone Jamming for Several 
Values of E(/Nq, Rb=10, Rj=10 

a consequence of the signal-to-jammer phase term in the bit error analysis. When the 
signal amplitude is consistently greater than the jammer’s amplitude, the phase difference 
has little effect and the rate of improving performance increases. This slope change is 
also noticeable in Figure 43 for the 20 dB curve. 
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Figure 43. Worst Case Multi-Tone Jamming for Several 
Values of Eb/No, Rb=100, Rj=l 

Figure 43 is a composite of performance under channel conditions favorable to the 
communicator. However, the 10:1 ratio in the signal's direct-to-diffiise over the jammer's 
direct-to-difflise ratio provides little relief to the communicator. 

Figure 44 displays the performance when the jammer's intelligence is relaxed to 
allow two interference tones in a hop slot when neither the signal nor the jammer 
experience fading. 
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Figure 44. Worst Case Multi-Tone Interference when 
Allowing Two Januning Tones per Hop Slot 

Figure 44 exhibits the same linear improvement in performance as 'E.JNj decreases 
for several values of Eb/No . There is a marked change in slope at ^JNj =0dB. Overall 
performance is better than the perffomance when at most one tone per hop slot is jammed 
in a near Gaussian channel as depicted in Figure 40. 
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IV. CONCLUSIONS 



A. FFH/MFSK PARTIAL-BAND NOISE JAMMING 

A performance analysis has been completed for the conventional fast 
frequency-hopped M-ary orthogonal frequency-shift keyed receiver employing diversity 
with noncoherent reception of signals transmitted over Ricean fading channels. We can 
draw several conclusions based on the results. 

Diversity alone is insufficient to overcome the effects of worst case partial band 
jamming acting in Ricean fading channels. The intelligent jammer can optimize the 
fraction of bandwidth jammed to force a worst case bit error ratio (BER). The 
noncoherent combining losses experienced by noncoherent systems is aggravated by the 
low signal-to-total noise ratio. These two factors combine to mitigate the advantages of 
increased diversity in a constant energy per hop system. 

Increasing the modulation order provides only modest improvements in performance 
which decrease as diversity increases. For most channel conditions, the performance for 
M=A and M=8 is essentially the same for all L greater than four. The only remaining 
advantage for increased modulation order is to offset the reduction in data rate that results 
from increasing diversity, Z,, in a constant energy per hop system. 

When the signal experiences Rayleigh or strong Ricean fading, the jammer need not 
be very sophisticated to cause marked impairment to the communicator. Additionally, 
when comparing the partial-band jamming strategy for M=2 and L=\ in which the signal 
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suffers moderate to strong fading to multi-tone jamming under the same conditions, we 
see that the partial-band jamming strategy causes a higher BER while requiring less 
precise knowledge of the receiver by the jammer. 

B. FFH/BFSK MULTI-TONE JAMMING 

In this thesis, the case of L=\ is considered for the fast frequency-hopped receiver 
described above. The influence of Ricean fading of both the signal and the jammer is 
included in the analysis as is the effect of thermal noise. 

First, the counter-intuitive result is obtained that Ricean fading of the jammer 
provides little, if any, relief to the communicator, while, as expected, moderate to strong 
fading of the signal tone markedly reduces performance. Worst case multi-tone jamming 
combined with signal fading produces unacceptably high probabilities of bit error, i.e., 
greater than 10'^. Above this level it is unlikely that forward error correction coding can 
be successfully implemented. 

The intelligent jammer can select the number of tones to jam which produces a worst 
case performance. The value of q chosen in the fading channel with thermal noise to 
cause worst case performance is the same value of q arrived at analytically for the noise 
free case without fading. For all combinations of channel conditions, not until the 
signal-to-jammer power ratio exceeds 1 OdB is the performance essentially thermal noise 
limited. When the jammer has a 3dB or better power advantage, the performance, in even 
the most optimistic channel conditions, is very poor and is dominated by the jammer 
power. 
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Relaxing the jammer intelligence to allow the case of possibly two interference tones 
in a hop slot is less efficient from the jammer's point of view. Performance is worse 
when the jammer is intelligent enough to jam not more than one signal tone per hop. 
Moreover, to be effective, a multi-tone jammer needs to be more intelligent than a partial 
band jammer since the former needs to know the exact location of the signal frequencies 
in addition to the bandwidth knowledge required of the latter. 

C. RECOMMENDATIONS 

The added complexity of a weighted diversity summer (noise normalized receiver) is 
warranted in most anti-jam systems. Efforts to reduce the uncoded BER below 10'^ will 
likely produce better overall performance when forward error correcting codes are 
implemented. When available, side information should be used to avoid a jammer or 
discard jammed hops. 

The partial-band noise jammer who is intelligent enough to force the worst case 
scenario is more threatening to the communication system than multi-tone interference 
and will require a detailed threat analysis for hostile sources. Though "dumb", 
non-hostile interference sources also warrant attention in mission planning since even a 
single signal tone consistently jammed can cause detriment to the communicator. 
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